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Abstract 

A beam splitter is an important component of an atomic/optical Mach-Zehnder interferometer. 
Here we study a Bose Einstein Condensate beam splitter, realized with a double well potential 
of tunable height. We analyze how the sensitivity of a Mach Zehnder interferometer is degraded 
by the non-linear particle-particle interaction during the splitting dynamics. We distinguish three 
regimes, Rabi, Josephson and Fock, and associate to them a different scaling of the phase sensitivity 
with the total number of particles. 
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I. INTRODUCTION 



Sub shot-noise interferometric measurements have become the subject of lively experi- 
mental and theoretical studies in view of possible breakthrough technological applications 
(for a recent review see In particular, interferometry with dilute Bose- Einstein conden- 
sates has become an important tool for experiments in fundamental and applied physics Q|. 



Among these, we mention the recent realization of double-slit 3^ Inland Michelson-Morley 
interferometers [5j], and the study of Hambury Brown-Twiss effect |6(. An archetypal two- 
mode interferometer is the Mach-Zehnder configuration, where two input fields are mixed 
in a beam splitter, undergo a relative phase shift 6, and are recombined in a second beam 
splitter. Two detectors, placed at the two output ports, allow the measurement of the total 
and relative number of particles. Single-atom detection with nearly unit quantum efficiencv 
has been recently demonstrated with Bose Einstein condensates in an optical box trap 3, |8( . 
From the collected data, it is possible to infer the value of 9 with a certain sensitivity which 



mostly depends on the nature of the input fields |9j. The goal of quantum interferometry 
is to detect a weak external phase shift with the maximum sensitivity. It has been shown 
[lol | that quantum mechanics imposes a fundamental uncertainty on the precision A8 with 
which the phase shift 9 can be measured. This ultimate limit of phase sensitivity is usu- 
ally discussed as the Heisenberg limit, A8 ~ being N? the total number of particles 
(atoms or photons) passing throu gh t he arms of the interferometer. Different schemes have 
been proposed to reach this limit I12I llH ll4|. In this report we focus on the Twin-Fock 
state first proposed in 



N N 



(1) 



This state provides the Heisenberg limit of phase sensitivity when it feeds the a and b input 
ports of a linear Mach-Zehnder interferometer. While it is very difficult to create the state 
(JTJ) with photons 



151 ] . Bose-Einstein condensates make possible the production of Twin- 



Fock states with a large number of particles through splitting an initial condensate using a 
ramping potential barrier. The transition from the superfluid to the Mott-insulator regime 
has been recently demonstrated in an array of wells [la, LLZ| • The dynamical splitting of a 
condensate into two parts has been experimentally studied in and theoretically analyzed 
in 19, 2(j] . Alternatively, the state in Eq.(^) can be created with two condensates realized 



independently 



2l|. 
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Here we analyze a crucial component of a Mach-Zehnder interferometer, namely the beam 
splitter which, in quantum interferometry, transforms an uncorrelated input quantum state 
to an highly entangled one necessary to overcome the shot noise limit 1 / tJNt [22, 23] . 
The beam splitter is created by a double well potential with a time-dependent barrier, 
taking into account the non-linear effects due to the particle-particle interaction in each 
condensate. Non linearity makes the condensate dynamics highly non-trivial. We show how 
non-linearity can degrade the sensitivity from the Heisenberg limit, in the non interacting 
case, toward the shot-noise limit, in the presence of strong interactions. We show that 
there is a range of values of the ratio between the non-linear interaction and the tunneling 
strength in which sub shot-noise sensitivity can still be achieved, provided the splitting is 
performed in the right time interval. The non-linear interaction can be tuned, for example, 
with a Feshbach resonance. A beam splitter for a Bose-Einstein condensate has recently been 
experimentally demonstrated, starting from a single condensate, using competing techniques: 
in atom chips with trap deformation fl Q, with matter wave Y-gnide Q and with a 
Bragg pulse [2a] . Among these, the trap deformation seems to be the most appropriate 
way to couple two independent condensates. Very stable optical double-well traps have 



recently been experimental 



y reported 0,li,|27|. Those have found applications in the study 



of Josephson dynamics [27| and matter wave coherent splitting ^, Q]. The developing of 
a beam splitter for Bose-Einstein condensates represents a challenging technological step 
toward the building of a matter wave ultrasensitive interferometer. 

In section II we will introduce our beam splitter model based on a two-mode approxi- 
mation of the double-well dynamical splitting. The relevant results of our analytical and 
numerical studies are presented in section IV. In particular we discuss how the phase sensitiv- 
ity, which, in the linear limit, is given by the Heisenberg limit, is degraded by the non-linear 
particle-particle interaction in each condensate. A detailed analysis allows us to distinguish 
three regimes, Rabi, Josephson and Fock. The transition between these is characterized in 
sections V and VI. 
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FIG. 1: Schematic representation of the beam splitter for a BEC in a double square well potential. 
The red line represents the potential barrier, the blue line is the condensate density distribution 
obtained by solving a one dimensional Gross-Pitaevskii equation. The beam splitter is realized in 
three steps: A) we start with two independent condensates, B) we allow a tunneling between the 
two wells for a certain time t by decreasing the height of the potential barrier separating the two 
condensates, and C) we rise the potential barrier to suppress the tunneling. 

II. BEAM SPLITTER MODEL 

In the linear limit of non-interacting bosons, the 50/50 beam splitter is represented by 
the unitary transformation |"L1 | 

\il> ut) = e- t *y-^—'\ip inp ), (2) 

a and b being annihilator operators for the two input ports; \ipin P ) and |Vw}> the input and 
output states of the beam splitter, respectively. A linear beam splitter for photons can be 
realized with a half transparent lossless mirror, while, in the case of ions it is given by a tt/2 
Raman pulse 

For interacting Bose-Einstein condensates in a double-well potential, the beam splitter 
transformation, in a two-mode model, can be written as 

|lM0> = eA^*'^xm«M*))*\^), (3) 



where 



H=-^^-a^aPb-K(t)(a r b + Pa) (4) 



is the two-mode Hamiltonian E c (t) being the charging energy, proportional to the 

particle-particle interaction in each condensate, and K(t) the coupling energy, representing 
the tunneling strength between the two condensates. The ratio between these two parameters 
can be controlled by tuning the interaction strength by a Feshbach resonance or, dynamically, 
by adjusting the height of potential barrier. In our model, a beam splitter for interacting 
Bose-Einstein condensates is created through a three stage process (see Fig. ((TJ)), A) we 
start from two independent condensates (E c (0) ^ 0, K(0) =0), as described by Eq. (JTJ) ; B) 
we allow a tunneling between the potential wells by decreasing the height of the potential 
barrier separating them {E c (t) ^ 0, K{t) ^ 0); and finally C) we suppress the tunneling by 
raising the potential barrier. We recover the 50/50 linear beam splitter Eq. (j2J) when the 
interaction is switched off, E c (t) = 0, and / +o ° dt K(t) = 7r/4 jsij. In figure (JTJ) we present 
a schematic representation of the beam splitter for BEC in a double square well potential. 
An optical box trap with single-atom detection capability has been recently experimentally 
realized in 3, • We point out that our analysis can be easily extended to condensates in a 
double well potential of arbitrary shape. 

We study the dynamics of a system of N particles by projecting its quantum state onto 
the Fock basis |n) = \N — n) a \n) b . In general, the output state |Vw) can be written as 

N 

|Vw(i)) = $> n (t)|n>, (5) 

n=0 

where the coefficients c n (t) are given by c n (t) = (n\ip ou t(t)) . If the condensate is initially in 
a Twin-Fock state, |^(0)) = \N/2)\N/2), then we have \c N/2 {0)\ 2 = 1, and |c n (0)| 2 = for 
n 7^ N/2 (see Fig.(J2jA)). To study the phase sensitivity we consider the operator 0] 

= ^10X01, (6) 

271 

where \<j>) are the normalized phase states \<p) — z^m=o e% ^ ^ 2+m -* \m) . The operator 

(P) has a positive spectrum and J Q 27r d0 E((p) = 1. Therefore, it defines a positive oper- 
ator value measure (POVM). For an arbitrary state \ipout{t)), the normalized probability 
distribution is 

P(<f>,t) = ^_ii(0|^(t))(^(t)|0) 



1 

2^ 



N 



,i (N/2+n) 4> 
n=0 



2 



(7) 
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FIG. 2: Relative number and phase distributions in the linear limit, plotted at different times: 
figures A and B, t = 0; figures C and D, t = ir/(SK); figures E and F, t = ir/(4K). Here N = 40, 
and K{t) is kept constant during the dynamics. 

n the linear case, Eq.((ZJ) coincides with the optimal quantum phase estimate proposed 
in An additional linear phase shift 6 simply displaces the whole phase distribution, 

thus providing a phase shift-independent probability distribution. In order to estimate the 
phase sensitivity, we calculate the distance between the two first minima on both sides of the 
central peak. This method, although used by other authors [3^, does not take into account 
the effect of the tails of the phase distribution, and it can only give qualitative results. 
The phase sensitivity has to be calculated with a rigorous Bayesian analysis of quantum 
inference as done, in the linear case, and for the whole Mach-Zehnder interferometer in 
where the effect of the tails of the distribution is discussed in detail. The initial Fock state 
has a flat probability phase distribution, corresponding to a complete undefined relative 
phase between the two condensates (see Fig. (013)). The linear beam splitter, for a constant 
K(t) = K, can be studied analytically, and the results do not depend on the particular value 
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of K ^ 0. The coefficients c n (t) are given by 



Cn(t) = e 



l N/2\N/2\ 
y (N-n)\n\ 



sin (tK)) 2 (cos(tK) 





[008(21*0], (8) 



where P" ,/3 [x] are the Jacobi polynomials. In general, during the dynamics, different c n are 
populated, and the phase distribution develops a central peak, as shown in the figures (J2JC) 
and (0D), referring to the linear evolution after a time t = n/(8K). In figures (J2JE) and 
(J2JF), we plot the \c n \ 2 and phase distribution after a tt/2 pulse (t = tt/(4K)) in the linear 
case. As we see, the spread of the \c n \ 2 distribution is of the order ~ N while the width of 
the main peak of the phase distribution is ~ 1/N. 

In contrast with the linear case, in non-linear dynamics, the time evolutions of K(t) and 
E c (t) play crucial roles. As a first approximation of the dynamical splitting, we consider a 
sudden displacement of the double-well at t = and for a certain time interval r. In this 
case, we have E c (t) = E c and K(t) = KQ(t — t), where Q(t) is the step function and K 
is a tunable parameter. We have checked, by 1-D Gross-Pitaevskii numerical simulations, 
that the parameter E c does not change significantly by displacing the double-wells, at least 
for a small overlap of the two-mode wave functions. This sudden displacement is the best 
scenario, as a slow separation of the wells will decrease the sensitivity, as shown in figure Q. 
In all our discussion, we develop a two-mode approximation, neglecting spurious excitations 
that can arise from the fast splitting of the wells. 

III. RESULTS 

With the formalism developed above, we now discuss the main results of our numerical 
and analytical study. First we analyze the dynamical change of the width of the phase 
distribution, keeping in mind that the smaller the width, the larger the phase sensitivity. 
The linear case can be studied analytically (see Eq. (|Hj)). and it is well known that the width 
attains its minimum at t — tt/(AK), corresponding to a tc/2 pulse or 50/50 beam splitter. 
The linear dynamics are periodic in time with period tc/(2K). In figure (@J) we present the 
width of the phase distribution as a function of time (in units tc/(AK)) for different values 
of the charging energy E c and for fixed values of K = 0.5 and N = 40. By increasing 
E c , the phase width reaches a minimum at smaller times and corresponding larger values. 
To these minimum of the phase distribution width there corresponds an optimal separation 
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FIG. 3: A) Dynamical evolution of K(t)/K(0) = (l - cos n (tvr/r)) , with r = 0.4 and for n = 2 
(yellow line), n = 4 (violet line), n = 6 (green line), n = 20 (blue line), n = 100 (red line) and 
n = oo (black line). B) Width of the phase distribution as a function of time for the different 
K(t) plotted in (A). The width of the phase distribution has a minimum at t ~ 0.4 for a sudden 
displacement of the double well (case n = oo, with K(t)/K(0) = 0(r — t)). For finite values of 
n (slow displacement) the minimum value of the width is reached at longer times but with larger 
values (lower sensitivities). In the insert we plot the phase distribution width as a function of time 
in logarithmic scale. The best scenario is obtained for the fastest splitting. Here we considered 
K = 0.5, AT = 40 and a constant E c = 0.4. Time is in units tt/(4K). 

time to realize the non-linear beam splitter. In Fig. (JH) we have plotted the dynamics just 
after the minimum. At longer times, the dynamics become almost chaotic, and the absolute 
minimum, together with the optimal separation time, is the one considered in the figure. 

In figure (0) we show the relative number and phase distribution corresponding to the 
maximum sensitivity for given E c . It can be compared with the linear case E c = (see 
Figs.(j2lE) and (0F))- As shown in Fig.fjSJ), for small values of E c the phase distribution 
matches the linear one (compare Fig.(0B) with Fig. (OF)), and it is characterized by a 
narrow central peak. We notice that the 7r-periodicity of the perfect linear case is lost 
in place of a 27r-periodicity (this effect characterizes the presence of a non linear coupling 
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FIG. 4: Width of the phase distribution Eg .([7]) as a function of time (in units 7r/(4X)) and for 
different value of the interaction energy E c . The linear case (E c = 0) is represented by the red line. 
By increasing E c the minimum is attained at smaller times and at corresponding larger widths. 
We notice that, independently of E c , at t = we have a flat probability distribution (see Fig.(|2jB)) 
corresponding to a phase width 2ir and complete phase uncertainty. Here K = 0.5 and N = 40. 

and will be discussed in the following). Increasing E c , the phase distribution broadens, 
as shown in Fig.fjHlD) and eventually becomes flat, as in Fig.(j3F). This matches a loss of 
phase sensitivity due to the non linear particle-particle interaction in each condensate. The 
|c n | 2 distribution becomes progressively narrow. When E c 3> K, for a general input state 
\i>mp) = Yln=o c n\n), we have 



The dynamical evolution of each c n is simply given by an evolution of its phase, corresponding 
to a time invariance of |c n (t)| 2 . The beam splitter becomes inefficient and it does not 
appreciably modify the input state. 

In Fig. © we present, for different values of E c and different numbers of particles, the 
optimal splitting time, which defines the splitting time giving the best phase sensitivity. In 
the linear limit it is given by t — tc/(AK) regardless of the numbers of particles. It decreases 
by increasing E c , and in the limit E c 3> K, when the dynamics are described by Eq.fJHI), the 
optimal splitting time becomes independent of N. This effect is highlighted in the figure by 



N 




(9) 



n=0 
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FIG. 5: Relative number and phase distributions at the optimal working point (minimum phase 
width), for different values of E c : A and B, E c = 0.04, C and D, E c = 4, E and F, E c = 400. These 
distributions can be compared with the ones in Figures (0E) and ifflF) calculated in the linear 
regime, E c = 0, and at the optimal time tt/{AK). Here the parameters are the same as in Fig.Q, 
N = 40 and K = 0.5. 

the asymptotic matching of different lines corresponding to different number of particles. 
Asymptotically in the number of particles, we can define the phase sensitivity as 

w = w- (10) 

where the prefactor a and the scaling factor (3 depend, in general, on the parameters K and 
E c . In figure (J7J) we show j3, as a function of K/E c and for different values of N (N = 40 in 
Fig.(J7jA) and = 80 in Fig.((7[B)). We can clearly distinguish between three regimes that 



32], 



characterize the two-mode dynamics. Following the notation introduced by Leggett 
we identify the Rabi, Josephson, and Fock regimes, depending on the ratio K/E c (see also 
). In order to find a qualitative definition of the three regimes, we consider the exact 



quantum phase model retrieved in 



By projecting the two-mode Hamiltonian (0J) 
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FIG. 6: Optimal splitting time (in units tt/(4K)) as a function of E c and for different values of N. 
Points are numerical results, lines are guides to the eye. Here K = 0.5. 

over the overcomplete Bargman basis we obtain the effective Hamiltonian 

E d 2 Kit) 2 
H eff (<p,t) = -y^r - K(t)N cos<P- cos 20. (11) 

The Rabi regime is defined by Kj E c ^> N. In this case the cos 20 term in the Hamiltonian 
(fTT]) is dominant over the cos (p. The effect of the cos 20 term is to dynamically squeeze 
the initial flat phase distribution creating a central peak with a width at the Heisenberg 
limit. In fact, in this regime, we have /3 ~ 1, corresponding to sub shot- noise sensitivity. 
In particular, for E c = 0, the cos0 term becomes negligible if compared with the cos 20 
term and the resulting phase distribution has period tt as noticed above (see Fig.(J2lF)). 
Physically, this is a consequence of the perfect symmetry of both the input Twin Fock state 
and the 50/50 beam splitter. As soon as E c ^ 0, the 7r-periodicity is lost as a consequence 
of the presence of the cos0 term in Eq. (jll|) . This effect can be observed in Fig.(0B). 
The Josephson regime is given by 1/N <C K/E c N. It is the dominant regime when 
we increase the number of particles, keeping fixed the ratio K/E c . As shown in figure (0), 
the phase sensitivity decreases when K/E c < N and, in the Josephson regime, we recover 
the shot-noise limit f3 = 1/2. The prefactor a becomes progressively large. In figures (0C) 
and (0D), we present the typical phase and \c n \ 2 distributions in the Josephson regime (in 
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FIG. 7: Scaling parameter (3 (denned in Eq. QlOj) ) as a function of K/E c . Points represent results of 
the numerical simulations, the blue line is a guide to the eye. The red vertical dashed lines delineate 
three regimes: Rabi, Josephson and Fock. Horizontal green dotted lines indicate the two relevant 
limits of quantum interferometry: the Heisenberg limit (3 = 1, and the Standard Quantum Limit 
[3 = 0.5. Note that the Heisenberg Limit is reached asymptotically in the number of particles. Here 
we considered A) N = 40 and B) N = 80. The Josephson regime becomes larger when increasing 
the number of particles. In the Fock regime it is not possible to define a scaling of the phase 
uncertainty since the phase distribution is almost flat and the phase sensitivity is A9 ~ 2ir. 

these figures, K/E c = 0.125, and N = 40): they have, to a good approximation, a Gaussian 
shape. 

The Fock regime is characterized by K/E c < and corresponds to strong interaction. 
In the case K/E c — > 0, the dynamics is described by Eq.<P|) and the phase distribution 
remains flat as shown in the Fig.(0F). It is not possible to define a scaling with N, and we 
have (3 — 0. In Figs. (0E) and (0F), we plot the narrow \c n \ 2 distribution and flat phase 
distribution characterizing the Fock regime (in the figures K/E c = 0.00125). 
In figure we highlighted the three regimes discussed above (dotted vertical lines). As 
the main result of this paper we see that, even in the presence of non-linearity, there is a 
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range of values of K/E c , corresponding to the Rabi regime, where it is possible to have sub 
shot-noise sensitivity, even in the presence of non-linear interactions. In the following two 
sections, we will discuss in detail the Rabi-Josephson and Josephson-Fock transitions. 

IV. JOSEPHSON-FOCK TRANSITION AND SELF TRAPPING 

We have observed a selftrapping effect associated with the two-mode dynamics of the 
system, which is different from the selftrapping of the Josephson oscillations It is 

possible to fully characterize this effect by examining the \c n \ 2 distribution. At t — the 
distribution is represented in Fig.(J2lA) where we have |ctv/ 2 (0)| 2 = 1. During the dynamics, 
I cn/2 (t) | 2 decreased, and the other modes |c n (t)| 2 with n ^ N/2 are populated. Since 
the energy is conserved, by monitoring the quantity |cjv/ 2 (t)| 2 , we can see how the initial 
energy dynamically distributes among all the modes. In particular, when |ca^/2 | 2 = 0, the 
initial energy has been completely distributed. In this way we can distinguish between two 
different behaviors: complete and incomplete energy distributions. When E c = 0, we have 
that \c]y/2(t)\ 2 follows a perfectly periodic motion with amplitude 1 and period r = ir/ (2K). 
For small values of E c , we find that, initially, the |cjv/2 (^) | 2 dynamics follows the linear 
behavior, and then it performs an almost chaotic motion. In this case the oscillations are 
still limited between 1 and ~ 0. For large values of E c we observe that |cAy 2 (t)| 2 oscillates 
between 1 and a value clearly different from 0, performing a selftrapped dynamics. The two 
different regimes, non- selftrapped and selftrapped are shown in figure (JBJ), where the time 
evolution of |c/v/ 2 (t)| 2 is reported. In Fig.Q we present the maximum oscillation amplitude 
of | Cat/2 {t) | 2 for different values of E c /K. In the non-selftrapped regime this amplitude is 
close to 1, while in the selftrapped region it is clearly smaller than 1. The selftrapped 
regime corresponds to an ineffective beam splitter where the two input modes are slightly 
mixed. This condition matches exactly the transition between the Josephson and the Fock 
regimes, where the beam splitter does not modify the input number and phase distributions. 
To find a qualitative estimate of the critical value (Ec/K)^f characterizing the transition 
between the two regimes, we make a three-mode approximation of the c n (t) dynamics. This 
corresponds to approximating Eq.© as \ip ou t{t)) = S^=/vy 2 _ 2 c n (i)|^), taking into account 
the symmetry of the distribution of c n around n = N/2 (See, for example Figs.(J2J) and (jHJ)). 
By substituting cjv/ 2 = ^fp^e 1 ^ , and analogously for ctv/ 2 +i and ctv/ 2+2 , in the limit N 3> 1, 
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FIG. 8: Time oscillation of |cjv/2 (^) 1 2 ? f° r different values of E c and in the selftrapping (A) and 
non-selftrapping (B) regimes. Time is in units tt/{AK). 

we obtain 

^c =2N y/poPlCOsOo - 0l) + y/pjpj COS(0 2 - 0i) 

With the three-mode constraint p^ + 2pi + 2p 2 = 1, where we consider p = 0.01, and after 
the direct numerical observation of p\ for the case N = 10, we obtain 

/ E \ JF 

-£) = 0.663 N. (13) 

V K / cr 

A simpler two-mode approximation would give (E c /K)^f = 0.899 N. The result is presented 
in Fig.© which shows a perfect agreement between the numerical results (points) and 
the three-mode approximation (line) given by Eq.([13p. We note that Eq.(|13|) marks the 
transition between the Josephson and the Fock regimes E c / K ^ N . 

V. RABI-JOSEPHSON TRANSITION 

To characterize the transition between the Heisenberg limit of phase sensitivity to the 
Standard Quantum limit, we introduce the entangled NOON state defined as 

_ |iV,0) + |0,iV) 

\*N00N) — ^ • 1-L4J 
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FIG. 9: Maximum oscillation amplitude of \cj^/2(t)\ 2 , as a function of E c /K and for different 
numbers of particles. We can clearly distinguish between the selftrapped and non-selftrapped 
regimes. Points are numerical results, lines are guides to the eye. Black points corresponds to 
N = 10, yellow N = 20, green N = 40, blue N = 80, red N = 160. In the insert we plot the value 
of (E c / K)^ corresponding to a maximum oscillation amplitude equal to 0.99, as a function of N. 
We define this as the critical point between the selftrapping and non-selftrapping regimes. Points 
represent numerical results with vertical bars giving the numerical uncertainty. The dotted line is 
the result of a three- mode approximation of the c n (t) dynamics, as given by Ea. (|13|) . 

When created after the first beam splitter of a Mach-Zehnder interferometer, this state 
leads to a Heisenberg limited phase sensitivity Qualitatively, this effect can be simply 
understood considering the probability distribution ((7j). For the NOON state (fl"4"|) c = cjv = 
1/V2 and we have P(<p) = cos 2 (iV0/2). This probability distribution is characterized by iV 
equal peaks in the interval [— 7r,7r], each with width 2ix/N. If the projection of the state 
IVw) (created from \ipinp) after the beam splitter) over the NOON state Eq. (fT4"|) is different 
from zero, then the state \ipi np ) can be used to reach the Heisenberg limit of phase sensitivity 
in a Mach-Zehnder interferometer |25j. From Eq.fjSJ), we can define the quantity 

CnqON = \(4>NOON\4>out)\ 2 = \co(t) +c N (t)\ 2 /2, (15) 
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FIG. 10: Cnoon-, defined in equation (|15|). as a function of E c /K. Points corresponds to results 
of numerical simulations using different number of particles: red points refer to N = 160, black 
N = 80, green N = 40, blue N = 20 and yellow N = 10. The lines are guides to the eye. The 
sharp decrease of Cnoon defines a critical value of E c /K which is reported, as a function of the 
number of particles in the inset. The dashed line in the inset is given by Eq. (|16[) and is obtained 
by numerical fitting. 

which gives the probability to obtain the NOON state, given the state after the first 

beam splitter. In general the quantity Cnoon depends on the parameters K and E c . For 
a Twin Fock input state we have co(i) = cjv(t) and Cnoon = 2|cAr(t)| 2 . The condition 
Cnoon - > marks the transition from the Heisenberg limit to the Standard Quantum 
Limit, thereby characterizing the Rabi-Josephson transition. In figure (fTUj) we plot the 
quantity Cnoon as a function of E c j K and for different numbers of particles. We notice a 
very fast decrease of Cnoon after a critical point (E c / K)f r J . In fig. (jl())l we plot this critical 
point as a function of the number of particles N. With a linear interpolation, we obtain the 
condition 

E \RJ 15 



KJcr N ' 

which exactly marks the Rabi-Josephson transition (K/E c ~ N) 
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VI. CONCLUSIONS 



Several aspects of interferometry with Bose-Einstein condensates have been discussed in 



the literature 



has been described in 



36l 137. 13811 - For instance, the effect of losses and a non-linear phase shift 



36| and |37[ respectively, and the detection efficiency, which seems 



to be the major obstacle toward the reach of the Heisenberg limit, has been studied in 
In this paper we analyzed how the non-linear effects associated with the particle-particle 
interaction in each condensate affect the realization of a BEC beam splitter. In particular 
we focused on two initially independent condensates in a Twin-Fock state. We showed, in a 
two-mode model, that the non-linear coupling decreases the interferometer phase sensitivity 
from the Heisenberg limit to the Standard Quantum Limit. Depending on the ratio Kj E c we 
characterized three regions for the phase sensitivity: the Rabi, Josephson and Fock regimes. 
We discussed the transitions between those regimes. The main result of our detailed analysis 
of the beam splitter is that there is an interval of the parameter Kj E c ^> N, the so called 
Rabi regime, where sub Shot Noise sensitivity can be achieved, despite the presence of a non 
linear coupling. This conclusion is of interest in view of recent experiments where both the 
particle-particle interaction (employing a Feshbach resonance) and the tunneling strength 
(tuning the potential barrier) can be appropriately controlled and changed, making the Rabi 
regime and sub-shot noise sensitivity achievable. 
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